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Abstract 

We present a new subtraction scheme for computing jet cross sections in 
electron-positron annihilation at next-to-leading order accuracy in perturbative 
QCD. The new scheme is motivated by problems emerging in extending the 
subtraction scheme to the next-to-next-to-leading order. The new scheme is 
tested by comparing predictions for three-jet event-shape distributions to those 
obtained by the standard program event. 



September 2006 



Contents 



1 Introduction 



2 Subtraction scheme at NLOl 



3 Notation 3 



3.1 Matrix elementsi 3 



3.2 Dimensional regularization. one-loop amplitudes and renormalizationi ... 4 



3.3 Cross sentionsl 6 



4 The approximate cross sectioni 6 



5 Integral of the approximate cross sectioni 9 



6 Checks! 14 



7 Conclusions! 14 



A Integrals over the factorised single-particle phase spaed 15 



A.l Collinear integrals! 15 



A. 2 Soft integraleJ 17 



B Calculation of symmetry factors! 18 



C Volume of the phase space in d dimensions! 20 



1 Introduction 



High-energy physics will enter a new era of discovery with the start of LHC operations in 
2007. The LHC is a proton-proton collider that will function at the highest energy ever 
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attained in the laboratory, and will probe a new realm of high-energy physics. The use of 
a high-energy hadron collider as a research tool makes substantial demands upon the the- 
oretical understanding and predictive power of QCD, the theory of the strong interactions 
within the Standard Model. 

At high any production rate can be expressed as a series expansion in as- Because 
QCD is asymptotically free, the simplest approximation is to evaluate any series expan- 
sion to leading order in as- However, for most processes a leading-order evaluation yields 
unreliable predictions. The next simplest approximation is the inclusion of radiative cor- 
rections at the next-to-leading order (NLO) accuracy, which usually warrants a satisfying 
assessment of the production rates. In the previous decade, a lot of effort was devoted to 
devise process-independent methods and compute rates to NLO accuracy and the problem 
is known to be solved PQ 12 El IH El E] ■ In particular, the dipole subtraction scheme 
provides a simple and fully universal way of computing the radiative corrections and has 
been implemented in two widely used programs the mcfm 13 |Hl E] and the nlojetH — h 
[ini CH 1121 codes. 

In some cases, typically and most importantly when the NLO corrections are large, 
the corrections at next-to-next-to-leading order (NNLO) accuracy are necessary in order to 
give a reliable prediction of the production rates. Recently, a lot of effort has been devoted 
to the extension of the subtraction schemes used at NLO to NNLO. It was found however, 
that the universal NLO schemes cannot be extended to NNLO I14j. which motivates 
the new method presented in this paper. 



2 Subtraction scheme at NLO 



The jet cross sections in perturbative QCD are represented by an expansion in the strong 
coupling as- At NLO accuracy we keep the two lowest-order terms, 

(y^-Lo = <y + o- . (2.1) 

Assuming an m-jet quantity, the leading-order contribution is the integral of the fully 
differential Born cross section dcrj^ of m final-state partons over the available m-parton 
phase space defined by the jet function J^, 

= I dalJm - (2.2) 

J m 

The NLO contribution is a sum of two terms, the real and virtual corrections, 

(T^^°= / da^+iJ„+i+ j dalJ^. (2.3) 
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Here the notation for the integrals indicates that the real correction involves m + 1 final- 
state partons, one of those being unresolved, while the virtual correction has m-parton 
kinematics, and the phase spaces are restricted by the corresponding jet functions J„ that 
define the physical quantity. 

In d = 4 dimensions the two contributions in Eq. ()2.3|) are separately divergent, but 
their sum is finite for infrared-safe observables order by order in the expansion in ag. The 
requirement of infrared-safety puts constraints on the analytic behaviour of the jet functions 
that were spelled out explicitly in Ref. jB]. 

The traditional approach to finding the finite corrections at NLO accuracy is to first 
continue all integrations to d = 4: — 2e {e ^ 0) dimensions, then regularise the real radiation 
contribution by subtracting a suitably defined approximate cross section dcr^'^^ such that 
(i) da^^i matches the point-wise singular behaviour of dcr^ in the one-parton IR regions of 
the phase space in any dimensions (ii) and it can be integrated over the one-parton phase 
space of the unresolved parton independently of the jet function, resulting in a Laurent 
expansion in e. After performing this integration, the approximate cross section can be 
combined with the virtual correction da^ before integration. We then write 



^NLO 



/ da^^^Jm+1 - da^^iJm + da^+ / da^'+i 



Jm . (2.4) 



J e=0 



R, A 

(Note that do"^'^;^ is multiplied by Jm, therefore, after integration over the phase space of 
the unresolved parton, it can be combined with do"^.) Since the first integral on the right 
hand side of Eq. ()2.4|) is finite m. d = A dimensions by construction, it follows from the 
Kinoshita-Lee-Nauenberg theorem that the combination of terms in the m-parton integral 
is finite as well, provided the jet function defines an infrared-safe observable. 

The final result is that one is able to rewrite the two NLO contributions in Eq. (|2.3|) as 
a sum of two finite integrals, 

-^^^ = / dar^« + / dar° , (2.5) 
that are integrable in four dimensions using standard numerical techniques. 



3 Notation 



3.1 Matrix elements 

We consider processes with coloured particles (partons) in the final states, while the initial- 
state particles are colourless (typically electron-positron annihilation into hadrons). Any 
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number of additional non-coloured final-state particles are allowed, too, but they will be 
suppressed in the notation. Resolved partons in the final state are labelled hj i, k,l, . . . , 
the unresolved one is denoted by r. 

We adopt the colour- and spin-state notation of Ref. [U]. In this notation the amplitude 
for a scattering process involving the final-state momenta {p}, abstract 
vector in colour and spin space, and its normalization is fixed such that the squared am- 
plitude summed over colours and spins is 

\Mm\' = {MJ\MJ . (3.1) 

This matrix element has the following formal loop-expansion: 

|A^) = |>l(o)) + |A^W) + ... , (3.2) 

where lA^^^-*) denotes the tree-level contribution, \Ai^^^) is the one-loop contribution and 
the dots stand for higher-loop contributions, which are not used in this paper. 

Colour interactions at the QCD vertices are represented by associating colour charges 
Ti with the emission of a gluon from each parton i. In the colour-state notation, each 
vector \A4) is a colour- singlet state, so colour conservation is simply 

(j2T?j\M) = 0, (3.3) 

j 

where the sum over j extends over all the external partons of the state vector and 
the equation is valid order by order in the loop expansion of Eq. ()3.2|) . 

Using the colour-state notation, we can write the two-parton colour-correlated squared 
tree amplitudes as 

l-^(°i)(M)r = {M^'\{P})\T,-T, \M^'\{p})) . (3.4) 

The colour-charge algebra for the product {Ti)"'{Tk)"' = Ti-Tj. is: 

Ti-Tk = Tu-Ti if i^k- T^i=Ci. (3.5) 

Here Ci is the quadratic Casimir operator in the representation of particle i and we have 
Cp = TR(Ar2 _ i)lN^ = (iV2 _ l)/(2iVe) in the fundamental and Ca = 2T^N^ = N, in the 
adjoint representation, i.e. we are using the customary normalization Tr = 1/2. 

3.2 Dimensional regularization, one-loop amplitudes and renor- 
malization 

We employ conventional dimensional regularization (CDR) in c? = 4 — 25 space-time dimen- 
sions to regulate both the IR and UV divergences, when quarks (spin-i Dirac fermions) 
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possess 2 spin polarizations, gluons have d — 2 helicity states and all particle momenta are 
taken as (i- dimensional. 



Turning to the renormalization of the amplitudes, let the perturbative expansion of the 
unrenormalised scattering amplitude 1^4^) in terms of the bare coupling gs = ^/^na^ be 



47r 



(3.6) 



where g is a non-negative integer and is the dimensional-regularization scale. The renor- 
malized amplitudes \M.rn) are obtained from the unrenormalized ones by expressing the 
bare coupling in terms of the running coupling as(/i^) evaluated at the arbitrary renormal- 
ization scale n\ as 



a>2^ = a,(/i^)/i^^57i 



1 - 



47r 



+ 0(«s(/^k)) 



(3.7) 



where (3o is the first coefficient of the QCD 13 function for nf number of light quark flavours. 



{3.t 



In Eq. ()3.7p . is the phase space factor due to the integral over the (c?— 3)-dimensional solid 
angle, which is included in the definition of the running coupling in the MS renormalization 
scheme,* 



Sr. 



(27r)'^-3 r(l - e) 



(3.9) 



We always consider the running coupling in the MS scheme defined with the inclusion of 
this phase space factor. 

The relations between the renormalized amplitudes of Eq. ()3.2|) and the unrenormalized 
ones are given as follows: 



5: 



"s(^r) /it e-l 
An 



q/2 



(3.10) 



An 



S7' 



q/2 



An 



R. 



-S;'l^f^'i\A^^^)-l^Sr\A^:^)^ . (3.11) 



After UV renormalization, the dependence on /i turns into a dependence on yU^, so the 
physical cross sections depend only on the renormalization scale /i^. To avoid a cumbersome 

*The MS renormalization scheme as often employed in the literature uses = {Airye"^'"^. It is not 
difficult to check that the two definitions lead to the same expressions in a computation at the NLO 
accuracy. 
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notation, we therefore set /x^ = in the rest of the paper. Furthermore, after the IR 
poles are canceled in an NLO computation we may set e = 0, therefore, the and 
factors that accompany the running coupling in the renormalized amplitude do not give 
any contribution, so we may perform the 

«s(^|)^?f ^_iy/' (^Al^s;'\ (<tM.y^'^^ (3.12) 



substitution in Eqs. (nrTTH) - (nrn| . 



3.3 Cross sections 



In our notation the real cross section dcrj^^_]^ is given by 

da^+i = AT J2 d0".+i(W) {MtU{p})\M^Zii{p})) , (3-13) 

where M includes all QCD-independent factors, X]{m+i} denotes a summation over all 
subprocesses, S^rn+i} is the Bose symmetry factor for identical particles in the final state 
and d0m+i({p}) is the d-dimensional phase space for m+1 outgoing particles with momenta 
{p} = {pi, . . . ,pm+i} and total momentum Q, 



d(l)m+l{pi, ■ ■ ■,Pm+l] Q) 



m+1 



n (^^+(^'^ 



(27r)"5W(pi + ---+p„+i-Q). (3.14) 



The virtual contribution daX is 



(3.15) 



In the rest of the paper, we define explicitly the approximate cross section dcr^',^^ and 
compute its integral (^cr^+l ■ 



4 The approximate cross section 

The construction of the suitable approximate cross section da^'^j^ is made possible by the 
universal soft and collinear factorization properties of QCD matrix elements ^3^]. In 
Ref. [T7] we introduced symbolic operators Cj^, that perform the action of taking the 
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collinear limit (pf or soft limit^ (p^ 0) of the squared matrix elements, respectively, 
keeping the leading singular term. Using this notation, we defined the formal expression 
A|A^^^j^P, that matches the singular behaviour of the squared matrix element in all the 
singly-unresolved regions of the phase space. 



^2 



\M^^UP^^Pr,...)\'. (4.1) 



This expression cannot yet serve as a subtraction term because it is defined precisely only 
in the strict collinear and/or soft limits. It has to be extended over the whole phase space, 
which requires an exact factorization of the m + 1 parton phase space into an m parton 
phase space times the phase space measure of the unresolved parton. 



d(f)m+li{p}) = d(f)mi{P}) [dpi] . 

With this phase-space factorization we define the approximate cross section as 

da^^, = d0„ [dpi] A\M^Zi\^ 
where ^ has the same structure as Eq. ()4.1|) . 



(4.2) 



(4.3) 



A\M^Zii{p})\' = Y. 



Y.l^^rm)+Ui{p})-J2^M{p}) 

We now define all terms on the right hand side of Eq. ()4.4j) precisely. 
The collinear counterterm C^j.{{p}) reads 

CMp}) = 8^c^sf^''-{MMpY'''^)\Pf.u^^^^^^ 



(4.4) 



(4.5) 



where the Pj,j^{zi^ri Zr,i, k±^i^r', £) kernels are the d-dimensional Altarelli-Parisi splitting func- 
tions as given in Ref. The momentum fractions Zi^r and 



and 



yiir)Q 

while the transverse momentum k\ ,■ ^ is 



VrQ 
y[ir)Q 



Ci,rPr Cr,iPi ~l" CirP- 



^1" 
ir 1 



Ci,r ^i,r 



<^iry{ir)Q 



Uir 



0-iry{iT)Q 



(4.6) 



(4.7) 



^For the precise definition of the coUinear and soft hmits, refer to Refs. |H1 117|. 

■I-Note in particular that the ordering of the flavour indices and arguments of the AltareUi-Parisi kernels 
has no meaning in the notation of Ref. |17j . 
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We used the abbreviations yir = SirjQ^ = 2pi ■ Pr/Q'^, y{ir)Q = ViQ + UtQ with yiQ = 
'^Vi ' Q / Q"^ ^ VrQ = 2pr- ' Q / and is the total four- momentum of the incoming electron 
and positron, while pf^ and air are defined below in Eqs. ()4.1()|1 and ()4.1H) respectively. This 
choice for the transverse momentum is exactly perpendicular to the parent momentum p^^ 
and ensures that in the coUinear limit the square of fc^j^ behaves as 



^± i r ~ SirZf^iZi^r , (4.8) 

as required (independently of Cir)- Choosing 

^± i r ^± i the collinear limit as can be shown by substituting the Sudakov parametriza- 
tion of the momenta into Eq. ()4.7|) (with properly chosen gauge vector). Note however, that 
in a NLO computation, fulfilling Eq. ()4.8p is sufficient to ensure the correct collinear be- 
haviour of the subtraction term and the longitudinal component that is proportional to 
Qr does not contribute due to gauge invariance of the matrix elements, so we may choose 
Cir = 0. The m momenta {p}^^'^^ = {pi, . . . ,pir, . . . ,pm+i} entering the matrix elements on 
the right hand side of Eq. ()4.5|) are 

p^ = Y^{p':+p^-(^^rQn, p':. = -r^p'^, n^z,r, (4.10) 

1 dir 1 O^ir 



where ^ 

Oiir = - y{ir)Q - \lylr)Q-'^yir ■ (4.11) 

The soft and soft-collinear counterterms S^{{p}) and C^^S^{{p}) are 

C^rSrMP}) = Svras/x^^-^T,^ IM^J^Kip}^'^^ ■ (4.13) 

Sir ^r,i 

If r is a quark or antiquark, S^{{p}) and Cj^^Sj.{{p}) are both zero. The eikonal factor in 
Eq. is 

SUr) = ^, (4.14) 

SirSrk 

and the momentum fractions entering Eq. ()4.13|) are given in Eq. ()4.6|) . The m momenta 
{p}^'^) = {pi, . . . ,pm+i} {Pr is absent) entering the matrix elements on the right hand sides 
of Eqs. dl^ and (ITOl) read 

K = Aj:[Q, {Q-Pr)/XrM/Xr) , U^T, (4.15) 
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where 

Ar=yr^^' (4-16) 

and ^ ^ ^ 

.,r — n , 2(K+Ky(K+K)^ 2K''K^ 

The matrix A^[K, K] generates a (proper) Lorentz transformation, provided = K ^ 
0. 

The momentum mappings introduced in Eqs. ()4.10p and ()4.15p both lead to exact phase 
space factorization in the form of Eq. ()4.2|) where [dpi] is a one-parton phase space times a 
Jacobian factor, 

[dpi;J = JjPr. {VYQ) -^^Mv'r) ■ (4.18) 

With our definitions for the momentum mappings the Jacobian factor depends on the 
number of hard final-state momenta. In the case of the collinear mapping in Eq. (j4.1Up . the 
Jacobian is 

Jt\Pr,P^r■, Q) = , — 6(1 - , (4.19) 

^ {y,rr + VTrQ - VrQY + 4l/^~(l - 2/-q) 

and that for the soft mapping of Eq. ()4.15|1 is 

J^\pr: Q) = A("-i)('^-^)-2 0(A,) . (4.20) 
In Eq. fl4.19p a^r is expressed in terms of the variable Pir, 

J {y-. + - VrQ? + ^y~r (l - VTrQ) - hjrTr +yTrQ- VrQ) 

= ^7^ ^ ■ (4.21) 

This concludes the definition of the approximate cross section in Eq. ()4.3p . Note that 
our do"^'"^ in general contains fewer subtraction terms than the dipole scheme. Furthermore, 
we can decrease the number of terms in Eq. ()4.4j) . because the symmetric treatment of the 
collinear subtractions is convenient for bookkeeping, but not essential in an actual computer 
code. 



5 Integral of the approximate cross section 

Next we evaluate the integral of the singly-unresolved approximate cross section over the 
one-parton unresolved phase space. Let us begin with integrating the collinear counterterm 
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Cj^. The transverse momentum as defined by Eq. ()4.7|1 is orthogonal to ptr, therefore, 
the spin correlations generally present in Eq. (j4.5|) vanish after azimuthal integration jH]. 
Thus when evaluating the integral of the subtraction term over the factorised one- 

parton phase space [dp^^I^{pr,pir', Q)], we may replace the Altarelli-Parisi splitting functions 

fifr their azimuthally averaged counterparts Pf-f^- Then 



[dptl{Pr,P^r■,Q)]Cd{p}) = ^ ^ J ^ ) C.^iVTr n', m- 1, s) TI \M^^\{p}(-^)\\ (5.1) 



a. 



2n " \Q 



/i 



2 \ e 



where y~Q = 2pir ■ Q and^ 



[5.2) 

The evaluation of these integrals is discussed in Appendix |X] and here we give only the final 
results. 



The function Cj^d/^^; n, e) depends on the momentum of the parent parton and the 
fiavours of the daughter partons. The three independent fiavour combinations are (we have 
'^iriVTrQ'^^^^) = Cri(z/^Q;n,e) and 0^-^(1/-^; n, e) = 0^3(1/-^; n, e)) 



Cgg{x;n,e) = x 



-2e 



Cqg(a;; n, e) — -^x 



x\e) 



l-e 



(5.3) 
(5.4) 



and 



Cgg(a;;n,£) = 2x 



-2s 



2 f l(-i) {x- e) - li°) {x- e)) + 1« (x; e) - 1^^) (x 



(5.5) 



The In \x; e) functions are given explicitly in Ref. [TH] . 

The expansion of Eqs. ()5.3|l - ()5.5p in powers of e is performed using the techniques of 
[m 1201 to obtain 



Cqg{x;n,e) 



-i + |---ln(x) + 0(£°) 



2e e 



C„Ax; n, e) 



Cgg{x;n,e) 



R 



Ca 



3e 



+ 0(£°) 



9 11 4 

- + - Hx) + 0{e') 

Se e 



(5.6) 
(5.7) 
(5.8) 



§We chose the dependence on m in the argument of this function as m — 1 because it is due to dependence 
on m in the Jacobian factor in Eq. H4.19|l . where it appears as m — 1 in the exponent. 
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The finite parts, not shown here, depend on n and can be easily found for any given n 
using the program of Ref. 1201 • We quote those for n = 2, 3 and 4 in Appendix 1X1 

Next consider the soft countertem. Defining 

£ (q^) ~ ^' ^) = -8™s/i^^ j Wlln{Pr] Q)] ]^Sik{r) , (5.9) 

we obtain 

Finally, integrating the soft-coUinear subtraction, Eq. ()4.13p we get 

/ [dpS:ibr;g)]QA(M) = ^S, (^)'cS(m- l,.)T,^|A<W({p}M)P, (5.11) 

with 

ZTT J J Sir ^r,i 

The evaluation of the integrals in Eqs. ()5.9|) and ()5.12j) is again discussed in Appendix El 
and here we give only the final results. 

The soft functions '^ikhJihyiQ^VkQ'^^^^) expressed with the standard beta and hy- 
pergeometric functions |21J, 



^.Avih VlQ. VlQ\ e) = -!^iL^iiLJ£) 5(1 _ 1 - e)Bil - 2e, n(l - e)) 

x-^2Fi (l,l,l-£,l-^^ 1 , (5.13) 

ViqUkQ \ ViqUkQ 

while 

' n{l - e){l - 2e) 



CS(n,e) 



£2 



B{1- e,l- e)B{l-2e,n{l- e)) . (5.14) 



Using the expansion 

z 2Fi(l, 1, 1 - e, 1 - 2) = [I + Li2(l - z) + 0{e^)\ (5.15) 
for the hypergeometric function, we find 

S..fe^, I/IQ, y-kQ-. ^> ^) = -4 - ; E ^ + ; In + 0(.°) , (5.16) 

^ fc=i VigykQ 

CS(n,.) = 1 + ^X^1 + 0(6°). (5.17) 
fc=i 
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Notice that 



^ikiyihyjQ.ykQ-^rn,e) + CBim,e) = - In ^^ + 0(5°). 

^ yiQVkQ 



(5.18) 



We quote the finite part of this expansion in Appendix El 



We are now in a position to calculate da^'^j^. Let us begin by recalling the form of the 
fully differential real cross section da^_^i given in Eq. ()3.13|) . Accordingly, the approximate 
cross section times the jet function is 



dcT^'+lJ„=Ar ^ d(j)m+l{{p})^ 



{m+1} 



X 



S{m+l} 



(5.19) 



In order to evaluate dcr^',^^ we first use the phase space factorization property of Eq. (j4.2p , 
then perform the integration to obtain 



1 a. 



/i 



J^da^^^Jm = AT ^ d0„({p}) ^7—7 7^Se(y^ 



{m+1} 



S{m+i} 27r \Q 



2 \ e 



r ij^r 



-C,^im-l,s)Ti\M^:^^{{P})\'U{p}) 



-(Sim-l,e)T',\M(^\{p})\'JU{p}) 



(5.20) 



This result is not yet in the form of an m-parton contribution times a factor. In order to 
rewrite Eq. ()5.20p in such a form we still need to perform the counting of symmetry factors 
for going from m partons to m + 1 partons, which is very similar to the counting in Ref. jH]. 
We give the details of the calculation in Appendix El Inserting equation Eq. ()B.6|) into 
Eq. we obtain 



,2 \ ^ 



dcrl''^t'Jm=^fJ2^M{p})^^sJ^^] JU{P}) 



X 



E 



C,(m-l,e)r?|7MS'({p})r + ES., 



[m 



m)\' 



, (5.21) 
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where we have introduced the flavour- dependent functions 

C, = C,,-CS, C, = ^C,, + nfC,,-CS. (5.22) 

We write the final result, dropping the m-parton jet function that appears on both sides 
of Eq. (ISI2H), as follows 

j^daX, = da^(^I{m-l,e), (5.23) 
where dcr^ is the Born-level cross section for the emission of m partons and 



(5.24) 



X 



E 



C^iViQ-, m - 1, e) T- + ^ ^ikiVik, ViQ, VkQ-, m-l,e) TiTk 

k^i 

is an operator acting on the colour space of m partons that depends on the colour charges 
and momenta of the m partons in |A^^P (different from the insertion operator of Ref. [H]). 
The notation on the right hand side of Eq. ()5.23|1 means that one has to write down the 
expression for dcr^ and then replace the Born level squared matrix element 

\MS = {Mt^Mt^) . (5.25) 

with 

{Mt^\I{m-l,e)\M^^^). (5.26) 

The parameter m — 1 in the argument of the insertion operator matches the arguments of 
the collinear and soft functions, not the number of coloured legs in the matrix element (in 
this case it is one less). Using colour conservation (T^ = —^k^iTiTk) to combine the 
Cj and Sj^ contributions, Eqs. ()5.(i|l - ()5.8j) and Eq. ()5.18j) . it is straightforward to check that 
our insertion operator differs from that defined in (7.26) of Ref. jH] only in finite terms, 

/(M; m -!,£)= /(M; e) + 0(£°) , (5.27) 

where 



- £ (^) ' E (r?^ + ^.^ E T.T.\ in 

^ ' i \ k^i / 



(5.2^ 



with the usual flavour constants 



7. = ^Cp , 7. = y ■ (5.29) 

It follows that /i^do"^'^^, as defined here, correctly cancels all e-poles of the virtual cross 
section dcr^. As a result, the integrand of the m-parton contribution, 

is finite and integrable in four dimensions (the potential kinematical singularities are 
screened by the jet function Jm)- This finite integrand is given in Appendix 1X1 
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6 Checks 



The cancellation of the singularities is a strong check on the correctness of the proposed 
scheme. We have performed such checks by approaching soft or collinear regions of the phase 
space from a randomly chosen point and computing the ratio of the (m + 1) parton squared 
matrix element and the subtraction terms. This ratio always approaches one. Further 
checks can be performed by comparing predictions for distributions to the predictions 
of other well-established computer codes for computing QCD jet cross sections at the 
NLO accuracy. Currently, our scheme is worked out only for colourless particles in the 
initial state, therefore, we decided to compare the predictions for the three-jet event-shape 
distributions thurst (T) and C-parameter in electron-positron annihilation, when the jet 
function is a functional 

J„(pi, ...,pn;0)= 5{0 - Osipi, . . .,pn)) , (6.1) 
with Osi^pi, . . . ,pn) being the value of either r = 1 — T or C for a given event (pi, . . . ,Pn)- 

Once the phase space integrations in Eq. ()2.5|1 are carried out, the NLO differential cross 
section for the three-jet observable O at a fixed scale Q takes the general form 

We computed the Bo(0) Born-level predictions as well as the Co{0) correction functions 
and found complete agreement with the corresponding tables of the benchmark calculation 
of Kunszt and Nason ||22| . 



7 Conclusions 

We have defined a new subtraction scheme for computing NLO corrections to QCD jet 
cross sections. For NLO computations the new scheme does not provide any particular 
advantage over the already existing methods and gives identical predictions. The need for 
the new scheme was motivated by studies in trying to extend the existing NLO subtraction 
schemes to computations at the NNLO accuracy. 

The new scheme is completely general in the sense that any number of massless final 
state coloured or colourless particles are allowed. It is worked out for processes without 
coloured partons in the initial state. The extension to deep-inelastic scattering and hadron 
collisions does not pose conceptual difficulties, but left for later work. 
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A Integrals over the factorised single-particle phase 
space 

In this appendix we compute the coUinear and soft functions Cj^, and CS and present 
their expansions relevant to three-, four- and five-jet production. 



A.l Collinear integrals 



We recall the definition of the collinear functions, Eq. ()5.2|) . from which one trivially gets 

where the azimuthally averaged splitting kernels are 

"1 - z 



+ 



2 



+ (1 - Zr)Zr 



l-e 

1 + (1 - Zr? 



1 Ziy j Ziy 



ezr 



(A.2) 
(A.3) 
(A.4) 



We parametrise the factorised one-particle phase space in terms of the variables air-, Uir 
and momentum fraction Zry Zry i « the latter being defined in Eq. ()4.(i|l . We find 



[dpS(PnP.; Q)] = (1 - a.rf-'^^'-^^~^^ {Q^Y~' ^ 



X 6 (yir - air {air + y^Q- "irZ/^g) j <^Oiir ^Vir d^r 

X {z+ - [y„ {z+ - Zr) {Zr - Z-)Y^ 

X 9(1 - air) ©(ttir) 0(^+ - Zr) Q{Zr - Z_) , (A. 5) 
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The limits of the Zr--integral are 

(+)/ ~ \ — 

''"^ a„{2a^r + yrrQ-airyrrQ) 

z^-\air,yir,y~Q) = ■ ^ r, (A.6) 

[air + y~Q - airy-QK^ocir + y~Q - airyrrq) 

with 2;*^+^ + z^^^ = 1. We now insert Eqs. ()A.2|) - ()A.5|) into Eq. ()A.1|) and obtain the results 
presented in Eqs. (j5.3|) - (|5.5|) . with 

p1 p1 pz+{a,y,x) 

''^'^ ll^\x] e) = da dy dz 5 {y - a (a + x - ax)) 

Jo Jo Jz^{a,y,x) 

X (1 - [z4a,y,x) - z^a^y, x)]-'+'^ 



X 



X 



x[y{z+{y,a,x) - z){z - z^{y,a,x))] ' -z . (A.7) 

These integrals are invariant under the z ^ 1 — z transformation, therefore, not all are 
independent; those with positive and odd n can be expressed with the others. For instance, 
ln\x;e) = ^ln\x; e). The integral over y is trivial by making use of the 6 function. 
After a complicated sequence of intergal transformations, the other two integrals can be 
transformed into known integrals. The results can be found in Ref. 



In an actual computation we need the expansion of the collinear functions in e to 0{e^). 
The pole terms are independent of m and are given in Eqs. ()5.6|1 - ()5.8|1 . Here we give the 
0(£°) terms, denoted by PinC^^. For n = m — 1 = 1 we need the functions only at a; = 1, 
where 

^nC^^(l; 1) = 7 - y , -HnC,,(l; 1) = ^ ^^""^3"^^ ■ (A.8) 

For n = m — 1 > 2 the results for arbitrary n are somewhat cumbersome combinations 
of elementary, 2F1 and 3F2 functions. Equivalent simpler expressions can be given in the 
following way: 

^2 



rinC^g{x] n) =cf + 2 (^In^ x + $(1 - x, 2, 2n - 1) + Li2(l - x) - — 

n—l ^ 

+ (^dnfl - 3^ Inx - ^ — 2^1 1,1 + i,l - X] 

1=1 ^ 



71—1 , 

+ yin T :2Fi{2n-l-i,l,2n-i,l-x), (A.9) 

i=0 

2T I 1 

nnC^^{x; = ^ I 4« + Inx - 2Fi{2n - 1, 1, 2n, 1 - x) 

2n — 1 X ( 2 — x\ . X 

+ -2Fi(^l,2n,2n+l,^) ). (A.IO) 

16 



Table 1: Constants used in the finite part of tlie expansion of the coUinear functions Cr, 



m n C^S cf dnfl dn,i dn,2 dn,3 



09 19 _11 _3 1 

^ 4 3 2 

4 o 89 _17 _8 3 1 

24 4 3 2 3 

c; 4 959 _277 _17 5 3 1 

"^ 360 60 5 3 5 6 



The finite parts of the colhnear functions for the gluon are not independent from the other 
two: 

J=inCgg{x; n) = 2J=inCgg{x; n) - —J=inCgg{x; n) ~ - . (A.ll) 

J-R <J 

The constants and dn,i can be found in Tabled The function s, a) is the Lerch 
transcendent [211, defined by the series 

oo ^ 

for which numerical codes for the evaluation exist While numerical codes for the 

hypergeometric functions also exist, it is actually faster to expand the 2F1 functions using 
their series expansion because for integer arguments the representation is a sum of polinoms 
and logarithms. 



A. 2 Soft integrals 

We recall the definition of the soft functions, Eqs. ()5.9j) and ()5.12|1 . from which one trivially 
gets 

S.k{y~f„mQ,y-,Q-,m- l,e) = {QT'^^ I [dpSi(p.;Q)] (A.13) 

i^e J yirykr 

and 

(S{m-l,e)=^^m-'^' !w%{Vrm-—- (A.14) 

Parametrizing the phase space with energy and angles, these integrals can be computed as 
done in Appendix B of Ref. [Hj (see also [21j), leading to Eqs. fj5.13|) and ()5.14p . Although, 
in Eq. (|5.22p we have combined the CS functions with the colhnear ones because these are 
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multiplied with the same colour factor, nevertheless one can always use colour conservation 
(T^ = — Xlfc^j TiTk) to combine the CS and S^^. contributions. Therefore, in presenting the 
0(e) terms in the e-expansion, we write only the finite term of their sum, which is simpler 
than the individual contributions, 

^n[S,,iyr„ y-Q, y.^, n) + CS(n)] = ^ + Li Jl - ^) + 2 In {-^\ E \ • (A.15) 

^ \ Vik J yyiQUkQ J 



In order to spell out the finite part of the m-parton contribution, da^^, we define the 
finite part of the one-loop amplitude as 

|A^^nM)) = -^^(W;^)l-M^^HW))+-Hn|A<«(M)). (A.16) 



Then 



+ ® /(m - 1, e)l_^ = A/" J] d0^ 1 2 Re{M^^\{p})\J^n\M^^\{p})) 

{m} ' \ 

H — -y 



J^JPm[S^fc(l/ifc,l/iQ,l/fcQ;m- 1) + CS(m- 1)] \M^^.^^.^f^){{p})\' 

k+i 



+Hn [C,(y,Q; m - 1) + CS(m - 1)] \Mt\{p})\' 



(A.17) 



B Calculation of symmetry factors 



Consider an m-parton configuration with mj quarks of fiavour f^fhf antiquarks of flavour 
/ and rrig gluons. From this configuration we can obtain an m -|- 1 parton configuration by 
changing 

(i) nig — > nig + 1 or (ii) m/ — > nif + 1, nif — > nif + 1, nig — > nig — 1 . (B.l) 
The ratios of symmetry factors corresponding to the two cases are 

^{m} _ ...nigl _ 1 



S{m+i} ...{nig + l)\ nig + 1' 



(B.2) 



_ ...nijlnijlniyl _ niy 



S{m+i} ...{nif + 1)!(% + l)!(m, - 1)! (m/ + 1)(% + 1) 
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We then have 



{m+l} i r^i {m} '^^™> ^ ^ \ i=qf r=g i=g r=qf 

+EE- + EE- + EE- 

*='?/ '"=9 «=9 ^=9/ «=9 »'=9 

+e""4 („,,;)^,,,/ ee-^ee--i-(b.4) 



Also 



{■m} 1 / \ J /\ J \i=qjr=qf i=qj r=qf 

EE--- = K + i) E EE--- = K + i) E 

^^...=K+i)j 5:5:...=K+i)j .... 

«=g/ r=g ir^q^ i=g r=qf j^^^^ 

EE--- = K + i)E---' 



^=9 r=g 



ir =g 

i=qfr=qf « 



EE-- ""^^',f'^" E-. (B.5) 

t=q;r=q; « 



thus we find 



E 5r;;:EE- = E'"sfT ^ -+ 5: ■ 



ir=qf,r=g ir=g,r=qf 



+^E ■■■+^E ■■■ + __E 

ir=qf,r=g ir=g,r=qf ir=g,r=g 

-E'"'^( E ...^ E ...). (B.a) 



ir=g,r=qf ir=g,r=qf 



The soft contribution to each sum in Eq. ()B.3|) is nonvanishing only if r is a gluon, so we 
have indicated the fiavour of r in each summation in Eq. ()B.6|1 . 
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C Volume of the phase space in d dimensions 



In this appendix, we present a simple derivation of the formula, obtained in Ref. for 
the volume of the phase space of m massless particles in d dimensions, which is a side 
product of the phase-space factorization presented in Eqs. ()4.18|) and ()4.2U|) . 

According to Eqs. (ji^ and (P?TH)) 

j d0„+i(g) = j d0^(Q)/(Q;m,£), (C.l) 
where, using Eqs. fj4.16|) and ()4.20|) . 

Working in the cm. frame, we parametrise the one-particle phase-space measure with the 
energy and angles. The integrand, (1 — yrq) = (1 — 2E/Q), depends only on the energy, 
therefore, 

HQ-,m,e) = ,^ l-E'-'dE 1-2^] , (C.3) 



(27r)^-iio 2 V Q 

where 

27r'^/2 

is the volume of the d-dimensional hypersurface. Introducing the new variable x = 2E/Q, 
the integral in Eq. ()C.3|) is readily obtained, 

r((rf-i)/2) 

We can get rid of the a/tt factors by using the identity. 



23-2d il-d)/2 

I{Q; m, e) = (Q')^ r((d-l)/2) ^ " 2' -^)(d- 2)]/2j . (C.5) 



Thus, we derived the following recursion relation: 

Y ( {m-l){d~^2) \ p /d-2\ 
d0™+l(Q) = 2-V/^(Q^)^ ^ ( lrn+l)!d-2) \ j ^'^-^^^ • ^^-^^ 

Starting from the known expression for the two-particle phase space 

0,(g) = 2W-^/^(g^)^i^ffl^, (C.8) 

and using the recursion relation in Eq. ()C.7|) . it is easy to obtain the result quoted in 
Ref. for = 4 - 2e. 
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